Mon. Not. R. Astron. Soc. 000, 000-000 (0000) 



Printed 1 February 2008 (MN KT E X style file vl.4) 



On the Motion of the Local Group and its Substructures 



S. Rauzy 1 and V.G. Gurzadyan 2 ' 3 

1 ANPE and Centre de Physique Theorique - C.N.R.S., Luminy Case 907, F-13288 Marseille Cedex 9, France. 

2 University of Sussex, Brighton BN1 9QH, UK. 

3 Department of Theoretical Physics, Yerevan Physics Institute, Yerevan 375036, Armenia (permanent address). 



qq > Received 

ON 



(N 



> 

O 

m 
m 
O 
oo 
On 

a: 

6 

C3 



13 



ABSTRACT 

The problem of the relative motion of the substructures of the Local Group of galax- 
ies revealed via S-tree method, as well as of the velocity of the Local Group itself, 
is considered. The existence of statistically significant bulk flow of the Milky Way 
subsystem is shown via 3D reconstruction procedure, which uses the information on 
the radial velocities of the galaxies, but not on their distances. Once the bulk motion 
of substructures is estimated, in combination with the observed CMB dipole we also 
consider the mean velocity of the Local Group itself. Assigning the Local Group the 
mean motion of its main substructures we evaluate its peculiar velocity in Milky Way 
frame V LG -*mw = (-7 ± 303, -15 ± 155, +177 ± 144) or 178 km s^ 1 toward galactic 
coordinates I = 245 and b = +85. Combined with CMB dipole Vmw-»cmb, we obtain 
Local Group velocity in CMB frame: Vlg^cmb = (-41 ± 303, -497 ± 155, 445 ± 144) 
or 668 km s _1 towards I = 265 and b = 42. This estimation is in good agreement, 
within 1 a level, with the estimation of Yahil et al (1977). 
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1 INTRODUCTION 



The peculiar motion of the Local Group (LG) seems to be- 
come one of the essential problems of observational cosmol- 
ogy, since it can provide at least crucial contraints on prop- 
erties of the local region of the Universe. In more optimistic 
evaluation, it can be the cornerstone for the linear theory 
of gravitational instabilities and, hence, for many aspects of 
the Big Bang cosmology, in general. 

The important question here is, obviously, the conver- 
gence of various dipoles, e.g. of optical galaxies and clusters 
of galaxies, IRAS galaxies, X-ray clusters, X-ray galaxies, 
etc., with the dipole of the Cosmic Microwave Background 
radiation, if the latter is caused by the Doppler effect. The 
reliable determination of the LG CMB dipole itself is the 
central point here. 

Though there is no absolute convergence, all the dipoles, 
in general, either are aligned in the direction of the defined 
CMB dipole or differ within limits, which presumably are 
affected by the choice of particular samples. For example, it 
is noticed, that the convergence level is improved, if deeper 
samples of extragalactic objects are involved, i.e. the matter 
distribution on more larger cosmic volumes is taken into 
account (Branchini&Plionis (1996)). 

While the result by Lauer&Postman (1994) on the 
quite differently oriented dipole obtained from the clus- 
ters of galaxies with velocities up to 15 000 km s _1 , 
seems to challenge the situation, the subsequent analysis 



by Branchini&Plionis (1996) involving even deeper survey 
of Abell/ACO clusters - up to 20 000 km s -1 , showed align- 
ment within 10° with respect to the CBM dipole, if the 
Virgocentric flow is also taken into account. 

The main common concern in these results based on the 
analysis of samples of various types of extragalactic objects, 
is in which degree the given population really traces the 
large-scale mass distribution. Though the general alignment 
of the dipoles of a given population with the CMB one, 
can imply the positive answer to this question, the reliably 
obtained divergence level of any of dipoles can be not less 
informative. Therefore, any alternative means of estimation 
of the motion of the LG should be of particular interest. 

In the present study we consider the possibility to es- 
timate the peculiar velocity of the LG using not the extra- 
galactic information, but its internal dynamical properties 
combined with data on CMB dipole. Namely, we analyse the 
substructure of the LG, and perform a procedure of 3D re- 
construction of the velocity of the whole system and its main 
subgroups. Thus, we continue to use the approach of 3D re- 
construction of the tangential velocities of galaxy systems, 
developed in Gurzadyan&Rauzy (1997), with the difference 
that in the present study the line-of-sight velocities have no 
components representing the Hubble flow. 

The first step is performed by the S-tree technique (see 
Gurzadyan&Kocharyan (1994)) enabling, in particular, the 
analysis of systems such as groups or clusters of galaxies. 
For the second step, we use the data on the distribution of 



© 0000 RAS 



2 S. Rauzy and V.G. Gurzadyan 



the line-of-sight velocities of the galaxies which are mem- 
bers of the LG and its subsystems, as revealed by the first 
step. Thus, we obtain the peculiar motions of the Milky 
Way and M31 subsystems with respect to each other. By 
the third step, from the bulk motions of the mentioned both 
subsystems we obtain the velocity vector of the Local Group 
in CMB frame. These results are discussed in term of the 
convergence of the various luminous dipoles with the one 
observed in the Cosmic Microwave Background Radiation. 



2 THE S-TREE TECHNIQUE 

The S-tree technique is developed for the investigation of 
properties of many-dimensional nonlinear systems and es- 
sentially uses the concepts of the theory of dynamical sys- 
tems. Refering for details to (Gurzadyan et al 1991, 1994; 
Gurzadyan&Kocharyan 1994 Bekarian&Melkonian 1997), 
here we outline its key points only. Its idea is based on the 
property of structural stability well known in theory of dy- 
namical systems enabling to study the robust properties of 
the systems with limited amount of information. The gravi- 
tating systems were known to be exponentially unstable and 
hence being among systems with strong statistical proper- 
ties. The advantage of the method in the context of galaxy 
clusters is in the self-consistent use of both kinematical and 
positional information on the clusters, as well as of the data 
on individual observable properties of galaxies - the magni- 
tudes. 

This approach is introducing the concept of the degree 
of boundness between the members of the given N-body sys- 
tem, i.e. the definition of a nonnegative function V called 
the boundness function which describes the degree of inter- 
action of a given subset Y of the initial set A with another 
subset X\Y according to some criterion. Particularly one 
can define a function 

V : S{A) -» R+ : (X,Y) ^ P X (Y), 

where S(A) is the set of all subsets of the initial set, when 
one of the sets contains the other one. The procedure of 
such splitting of A can be measured by a non-negative real 
number p using the boundness function, so that the problem 
is reduced to the finding out of a function E(p) denoting the 
set of all possible p-subsystems {.4i, . . . , Ad} of (A, V). 

Attributing for the given p the matrix D to another 
matrix T in a following way: 

T ab = if D ab < p, D ba < p, . . 

r ab = 1 if D ab > p, D ba > p; 1 > 

the problem of the search of a p-bound cluster is reduced 
to that of the connected parts of the graph F(p) - a S-tree 
diagram. 

For the matrix D ab a representation via the two- 
dimensional curvature K of the phase space of the system 
can be used: 

D ab =max{-^,0}, 

where p = (a, i),v = (b,j). The two-dimensional curvature 
is represented by the Riemann tensor 7? via the expression: 
K^ v — R^ v \ p u x u p , where u v is the velocity of the geodesies; 
the explicit expression of two-dimensional curvature for N- 
body gravitating system is derived in (Gurzadyan &Savvidy, 



1986) and has rather complicated form to be represented 
here. The advantage of the use of geometric characteristics 
such as Riemann, Ricci tensors, is well known in theory in 
dynamical systems, and has been used in astrophysical prob- 
lems as well. 

As a result, the degree of boundness between the mem- 
bers and subgroups of the given N-body system can be 
obtained, thus revealing the physically interacting system 
(cluster) and its hierarchical substructure. 

The computer code based on the S-tree method has 
been used for the study of the substructure of the Local 
Group of galaxies (Gurzadyan et al. (1993)), of the core of 
Virgo (Petrosian et al. 1997) and Abell clusters (Gurzadyan 
& Mazure, 1997). In these studies the information on the 
masses of galaxies has been used via the mass-to-luminosity 
ratio M — const L n , n = 1, though other relations - (n = 
0, 1/2), have been checked as well, and the robust character 
of the results on the subgrouping has been revealed; note, 
that these assumptions take into account also the existence 
of the dark matter associated with the galaxies. 



3 BULK FLOW RECONSTRUCTION 

The reconstruction of the 3D velocity distribution function 
from its observed line-of-sight velocity distribution is an in- 
teresting inverse problem. It was analytically solved by Am- 
bartsumian (1936) for stellar systems without any a priori 
assumption on the form of the soght function. The only as- 
sumption made was the independence of the distribution 
function on the spatial regions (directions), or equivalently, 
that the 3D velocity distribution <f>(vi,V2,V3) dv\ dv2 dv3 of 
such systems was invariant under spatial translations. It 
means, that the theoretical probability density dP t h of the 
system reads as follows 

dPth = <t>(vi, V2, V3) dvi dv2 dv3 p(r,l,b)r 2 cosbdr dldb (2) 

where <^>(t>i, V2, V3) is the 3D velocity distribution function 
(in galactic cartesian coordinates) and p(r,l,b) is the 3D 
spatial distribution function (in galactic coordinates). Un- 
der this assumption, Ambartsumian had proven the theo- 
retical possibility of reconstructing of 0(«i , «2, V3) from the 
observed radial velocity distribution function f(v r ,l,b) 

dPobs = f{v r , I, b) cosbdldbdv r (3) 

where v r = v.f is the projection of the velocity v = 
(vi,V2,V3) on the line-of-sight direction f = (fi, r~2, ^3) = 
(cos I cos b, sin I cos 6, sin b) . 

However, computer experiments show that the direct 
application of the Ambartsumian's formula is hardly possi- 
ble for that purpose. It is natural, since the derivation of a 
smooth function based on discrete information on relatively 
small number of particles (say, less than 1000) is a nonlinear 
problem. The number of points should exceed essentially the 
number of galaxies of real clusters in order to apply this for- 
mula with success. This fact is a consequence of the principal 
difference between the N-body problems in stellar dynamics 
and dynamics of clusters of galaxies. 

Nevertheless, some quantities of interest such as the first 
order moment of this distribution, i.e. the mean 3D veloc- 
ity of the system, can be obtained. Additional hypotheses 
on the distribution function 0(«i , «2, V3) are then necessary. 
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Hereafter, we assume that (p(vi,V2,V3) can be written as 
follows 

3 

4>{vi,v 2 ,v 3 ) dvidv 2 dv 3 = Y[ <t>i(yi) dv t (4) 

i=l 

where the distribution functions <j>i are centered on Vi and 
of variances Ui (4>i(vi) = (f>i{vi] Vi, 0»)). The statistics of 
the mean 3D velocity, i.e. the bulk flow velocity, are de- 
rived in the appendix A. The use of the maximum liklehood 
technique forces us to entirely specify the functions fa. We 
choose the Gaussian representation (fa = fa — fas = g 
with g Gaussian) and assume the isotropic velocity diper- 
sions (cri = (72 = (73 = <r„). The velocity field of the system 
is, thus, split into a mean 3D velocity vb = {v x ,v y ,v z ) = 
(vi,V2,V3) (i.e. a bulk flow) plus a 3D random component, 
isotropic Maxwellian, centered on and of velocity disper- 
sion a v . 

The fact that gravitating N-body systems do posess 
strong statistical properties peculiar to Kolmogorov systems 
(Gurzadyan&Savvidy, 1986; see also Gurzadyan&Pfenniger, 
1994), can be the justification for the separable form of dis- 
tribution functions in Eqs.(2) and (4). Note, that the rep- 
resentation (4) is valid if the potential is time independent, 
and the energy is an isolating integral. The exact Maxwellian 
form of the distribution, however, is not required for the 
analysis, its representation in the form g(vi,Vi,a) is suffi- 
cient. 

The derived bulk flow statistics is a robust estimator. 
It does not depend on the distances of galaxies under con- 
sideration. This is a positive feature, since there is still dis- 
agreement between various authors in the distance estimates 
of the members of LG. Moreover, it is shown in appendix 
A that the v x , v y and v z estimators do not depend on the 
projected angular distribution function r)(l, b) of the sample, 
even if their accuracies do. 



4 APPLICATION TO THE LOCAL GROUP 

4.1 The Local Group subsystems 

In the S-tree study of the substructure of the LG 
(Gurzadyan et al. (1993)) a sample of 39 galaxies has been 
used, compiled from various studies. The membership of 
galaxies in two main subgroups - dominated by the Milky 
Way and of M31, has been revealed, in general, confirming 
the conventional views. Physical connections between some 
individual galaxies have been also indicated, not reported 
before (e.g. of NGC 6822 and IC 1613). Some galaxies ap- 
peared to have no actual influence on the dynamics of LG 
and vice versa, and therefore, were considered not to be 
its members. The LG as a physical system of the influenced 
galaxies only, was shown to extend on less than 2 Mpc. These 
galaxies and their membership are listed Table 1. The recent 
studies (see Karachentsev (1995) and references therein), 
though record some differences in the data on the galaxies, 
in general, no drastic reliable change is observed. Existence 
of galaxies obscured by the Galactic disk, like newly discov- 
ered Dwingeloo 1,2 ones (Burton et al. (1996)), also cannot 
be ruled out, though it seems unlikely that, at least the lat- 
ter ones can essentially influence the substructure of the LG 
(Lynden-Bell, 1996). 



Table 1. Local Group data used for this analysis: the equato- 
rial coordinates, radial velocity v r (WITH RESPECT TO THE 
MILKY WAY REST FRAME) and their distance r, are listed for 
32 objects within 2 Mpc. M31 and Milky Way subsystems were 
deduced by using the S-tree method (Gurzadyan et al. (1993). 
Some typographical mistakes occured in Table 1 of the mentioned 
paper are corrected in the present Table). 



LOCAL GROUP : 32 galaxies 



Object a S v r r 

(1950) (1950) (km s" 1 ) (Mpc) 



Milky Way subsystem : 12 galaxies 



Milky Way 


00 00.00 


00 00.0 





0.0 


SMC 


00 51.00 


-73 06.0 


91-1-^ 


U.UDo 


Sculptor S 


00 57.60 


-33 58.0 


74±? 


0.085 


IC1613 


01 02.22 


01 51.0 


-153±3 


0.64 


Fornax D 


02 37.84 


-34 44.4 


-51±? 


0.19 


LMC 


05 24.00 


-69 48.0 


130±? 


0.052 


Carina 


06 40.40 


-50 55.0 


14±? 


0.093 


Leo I 


10 05.77 


12 33.2 


177±? 


0.22 


Leo II 


11 10.83 


22 26.1 


15±? 


0.22 


UM 


15 08.20 


67 18.0 


-88±? 


0.065 


Draco 


17 19.40 


57 58.0 


-95±? 


0.075 


N6822 


19 42.12 


-14 55.7 


55±5 


0.62 




M31 subsystem : 7 gal 


axies 




N147 


00 30.46 


48 13.8 


36±50 


0.7 


N185 


00 36.19 


48 03.7 


67±22 


0.72 


N205 


00 37.64 


41 24.9 


49±11 


0.64 


M31 


00 40.00 


40 59.7 


-110±1 


0.67 


M32 


00 39.97 


40 35.5 


86±6 


0.66 


LGS-3 


01 01.20 


21 37.0 


-100±8 


0.72 


M33 


01 31.05 


30 23.9 


-39±0.5 


0.82 


Others : 13 gal; 


ixies (with r 


< 2 Mpc) 




IC10 


00 17.69 


59 00.9 


-144±2 


1.25 


N55 


00 12.40 


-39 28.0 


97±3 


1.4 


IC342 


03 41.95 


67 56.4 


221±? 


1.84 


Leo A 


09 56.53 


30 59.2 


-17±10 


1.59 


N3109 


10 00.80 


-25 54.8 


193±? 


1.6 


Sex A 


10 08.57 


-04 27.7 


154±5 


1.3 


DD0155 


12 56.20 


14 29.2 


183±5 


1.0 


Sagittarius 


19 27.01 


-17 47.0 


10.9±1.4 


1.11 


DDO210 


20 44.13 


-13 02.0 


-13±10 


1.5 


IC5152 


21 59.60 


-51 32.0 


78±15 


1.51 


Hog 


23 23.80 


-32 40.0 


75±5 


1.3 


Pegasus 


23 26.05 


14 28.3 


-14±7 


1.0 


W-L-M 


23 59.40 


-15 44.6 


-56±8 


1.0 
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Table 2. Estimates of the bulk flow and of the velocity disper- 
sion for the Milky Way subsystem. Residuals are listed for each 
object. Standard deviations are computed by using numerical sim- 
ulations. 



Milky Way subsystem : 12 galaxies 



Bulk flow : vg = (v x ,v y ,v z ) 



Galactic coordinates Equatorial 
(cartesian) (1950) 

v x = 104 ±64 km s -1 Z = 312 a = 13 38.75 
Vy = -116 ±42 km s" 1 6 = 16 8 = -45 43.8 

v z = 45 ±33 km s" 1 ||v B || = 162 km s" 1 

Velocity dispersion : a v = 69 ± 15 km s — 1 



Table 3. Estimates of the bulk flow and of the velocity dispersion 
for the M31 subsystem. Residuals are listed for each object. Stan- 
dard deviations are computed by using numerical simulations. 



M31 subsystem 


: 7 galaxies 


Bulk flow : vg 


= (Vx,Vy,V z ) 


Galactic coordinates 


Equatorial 


(cartesian) 


(1950) 


v x = -117 ± 541 km s" 1 I 


= 144 a = 11 54.35 


v y = 85 ±267 km s" 1 b 


= 65 8 = 50 07.09 


v z = 309 ± 254 km s" 1 


|]v s || = 341 km s" 1 



Velocity dispersion : a v = 78 ± 20 km s 1 



Residuals 

Object v r v r — r.vg 

cr* =58 km s _1 



Residuals 

Object V r V r — r.Vg 



a* =57 km s 1 

N147 36 -16 

N185 67 15 

N205 49 36 

M31 -110 -120 

M32 86 78 

LGS-3 -100 -2 

M33 -39 



SMC 


-21 


-99 


Sculptor S 


74 


107 


IC1613 


-153 


-37 


Fornax D 


-51 


-23 


LMC 


130 


45 


Carina 


14 


-55 


Leo I 


177 


76 


Leo II 


15 


-25 


UM 


-88 


-21 


Draco 


-95 


-31 


N6822 


55 


31 



The choice of membership criteria is always crucial 
while studying the substructures of clusters of galaxies. In 
Karachentsev (1995), only closed Keplerian orbits are con- 
sidered as criteria of a membership, while strictly speaking, 
precise Keplerian closed orbits are never typical in few-body 
problem, as one has in the case of LG. As it is mentioned 
above, the S-tree technique enables to reveal the degree of 
influence, whatever the exact shape of the orbits is. 

Our aim here is to show the principal possibility of the 
reconstruction of 3D vector of the LG based on the substruc- 
ture analysis, having in mind, that the latter can be easily 
revised, if some observational input data on galaxies should 
be essentially modified with high enough confidence level. 

4.2 Bulk flow estimate 

The MW and M31 subsystems contain respectively 12 and 7 
galaxies. Such a small number of objects forbids indeed the 
use of the Ambarsumian's formulae in order to reconstruct 
the 3D velocity distribution functions of these two systems. 
Neverthless, in this section we show that it is possible to 
evaluate their 3D mean velocity, i.e. their bulk flow. 

For each subsystem, as it is assumed Eq. (2), the ve- 
locity distribution function <f)(vi, Wa, Vz) is translation in- 
variant, so that (j)(v\,V2,vz) can be split into a bulk flow 
vb — (v x ,v y ,v z ) plus a 3D random component, isotropic 



Maxwellian, centered on and of velocity dispersion a v (see 
section [| Eq.(4)). According to the bulk flow statistics de- 
rived in appendix A, Tables 2 and 3 give the values of the 
bulk flow estimates for the MW and M31 subsystems, re- 
spectively. These velocities are expressed in km s _1 and with 
respect to rest frame of the Milky Way (MWR frame). 

For each system, the values of the velocity dipersion a v 
and of the accuracy of the bulk flow estimate are calculated 
by using numerical simulations. Large number of samples 
have been simulated according to the following characteris- 
tics: 

- The angular position of each object of the simulated 
system is identical with the observed one. 

- The 3D velocity of each object is the estimated bulk 
flow (for example, v_b = (104,-116,45) in km s for the 
MW subsystem) plus a 3D random component, Maxwellian 
and isotropic, of velocity dispersion a v . 

- The simulated radial velocity of each object is the 
line-of-sight projection of its 3D velocity, plus a white noise 
accounting for measurment errors on v r (see table 1). When 
this information is missing, a value of 10 km s _1 is adopted. 

Since the quantity a v is unknown, we adopt the value 
of the biased estimator a* of a v as a fiducial starting point 
(see appendix A). After few iterations, the correct value of 
the velocity dispersion a v is reached when the a% estimate 
of the simulated samples corresponds to the observed one. 
Errors bars on the v x , v y , v z and a v estimates are afterwards 
calculated by using a large number of simulated samples. 

The effect of the Hubble expansion have been also in- 
vestigated. The radial velocities listed in table 1 are indeed 
not corrected for the Hubble flow. We have applied this 
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Table 4. Estimates of the bulk flow and of the velocity disper- 
sion for the Local Group, ignoring its substructures. Residuals arc 
listed for each object and the cr* estimates are given for the M31 
and Milky Way subsystems and for the 13 remaining galaxies. 
Standard deviations are computed by using numerical simula- 
tions. 



LOCAL GROUP : 31 galaxies (with r < 2 Mpc) 



Bulk flow : vj = (v x ,v y ,v z ) 



Galactic coordinates Equatorial 
(cartesian) (1950) 



v x = -28 ± 45 km s 


,-i 


1 = 233 a = 08 13.35 


Vy = -37 ± 32 km s 


r 1 . 


b = 13 5 = -11 00.53 


v z = -11 ±30 km e 


,-i 


vs|| = 47 km s _1 


Velocity dispersion : 


: (Jy — 103 i 14 km s 1 




Residuals 


Object 




V r V r — r.VB 


Milky Way subsyst 


cm : a* = 73 km s _1 


SMC 




-21 -25 


Sculptor S 




74 86 


IC1613 




-153 -138 


Fornax D 




-51 -60 


LMC 




130 109 


Carina 




14 -21 


Leo I 




177 78 


Leo II 




15 13 


UM 




-88 -76 


Draco 




-95 -69 


N6822 




55 98 


M31 subsystem 


: cr* = 78 km s _1 


N147 




36 56 


N185 




67 87 


N205 




49 69 


M31 




-110 -90 


M32 




86 106 


LGS-3 




-100 -83 


M33 




-39 -27 


Others (with r 


< 2 Mpc) : cr* = 121 km s" 1 


IC10 




144 -163 


N55 




97 110 


IC342 




221 222 


Leo A 




-17 -49 


N3109 




193 151 


Sex A 




154 113 


DD0155 




183 170 


Sagittarius 




10.9 52 


DDO210 




-13 30 


IC5152 




78 98 


Hog 




75 97 


Pegasus 




-14 19 


W-L-M 




-56 97 



correction by substracting to the radial velocity v r of each 
objects its Hubble velocity Hor. Since the LG galaxies are 
close to the Milky Way, the bulk flow estimate remains al- 
most unchanged. For H = 100 km s _1 Mpc -1 , the MW 
bulk flow becomes vj = (106,-128,60) (i.e. vb = 176 km 
s _1 toward I = 310 and b = 20) and the M31 bulk flow 
■v B = (-2, 77, 327) (i.e. v B = 336 km s _1 toward I = 91 and 
6= 77). 

In Table 3, bulk flow components and associated un- 
certainties of the M31 subsystem are expressed in galactic 
cartesian coordinates. In the coordinates frame such that X- 
axis is aligned with M31 line-of-sight (i.e. toward I = 121.2 
and 6 = —21.6) and Y-axis lies in the galactic plane, the 
3D bulk flow rewrites v s = (10 ± 34, 72 ± 623, 335 ± 291). 
In this coordinates, bulk flow of the MW subsystem reads 
v s = (-155 ± 57, -26 ± 53, -15 ± 26). 

In Table 2, one can notice that 3D bulk flow estimates 
of the MW subsystem does not take into account the ra- 
dial velocity of Milky Way itself. Since the analysis is per- 
formed in the MW rest frame, the MW velocity is zero by 
definition. However, no prefered line-of-sight direction can 
be assigned to the Milky Way by a MW observer. In order 
to overcome this principal difficulty, let us proceed to the 
following thought experiment. The observer ventures a foot- 
step outside Milky Way, small enough such that the radial 
velocities of others galaxies remain unchanged. It allows in- 
deed to define a line-of-sight direction for MW, opposite to 
the footstep walk. The bulk flow statistic derived appendix 
A is then applied to the MW subsystem, including the zero 
radial velocity of MW with its associated direction. This 
scheme is repeated for a large number of random footstep 
directions, such that the MW line-of-sight direction distri- 
bution becomes isotropic. Averaging over all these bulk flow 
estimates finally gives the 3D mean velocity of MW subsys- 
tem. The results of Table 2 remain almost unchanged, i.e. 
v s = (81±68, -104±44, 38±34) or v B = 137 km s" 1 point- 
ing toward I = 308 and b — 16, with a velocity dispersion of 
71 ± 16 km s" 1 . 

We have also performed our analysis on the 32 Local 
Group galaxies situated nearer than 2 Mpc, ignoring the 
MW and M31 dynamical substructures. Bulk flow and veloc- 
ity dispersion estimates are shown in Table 4. Note, however, 
that these estimates have no much sense since the presence 
of the MW and M31 subsystems rules out our main assump- 
tion, i.e. the velocity distribution of the 32 Local Group 
galaxies is not invariant under spatial translations. The val- 
ues of the biased velocity dispersion estimator a* for the 
MW and M31 subsamples are more interesting. They are 
significantely greater than the a* estimates of tables 2 and 
3. This fact indicates the existence of kinematic substruc- 
tures in the Local Group, and so, strengthens the present 
analysis. 



5 CONCLUSION AND DISCUSSION 

The motions of the two subsystems of the Local Group have 
been estimated. The M31 and MW dynamical subsystems, 
containing respectively 7 and 12 galaxies, have been iden- 
tified via the S-tree technique, which takes into account 
in a self-consistent way the degree of influence of one ob- 
ject (or a set of objects) on another, whatever their or- 
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bits can be. In the rest frame of the Milky Way (MWR 
frame) the bulk flow statistic derived appendix A, gives 
v B = (104±64, -116±42, 45±33) in km s _1 and in cartesian 
galactic coordinates (or vb = 162 km s _1 pointing toward 
I — 312 and b = 16) for the motion of the MW subsystem, 
and v B = (-117±541,85±267, 309 ±254) or v B = 341 km 
s _1 toward I — 144 and b = 65 for the motion of the sub- 
system dynamically associated with M31. While these esti- 
mates have been derived from the radial velocity and angu- 
lar position of galaxies, the information on the distance has 
not been used (the point has its importance since the main 
ambiguity in the data concerns the galaxy distances). Note, 
that the independence on the distances concerns only the 
velocity reconstruction procedure, while the information on 
distances is used for S-tree analysis; however, the statistical 
results of the latter concerning the subgrouping properties 
are robust relative the error-boxes of data, unless some data 
will be modified drastically. Note the following two points. 
On one hand, the relative velocity of M31 subsystem with re- 
spect to the MW subsystem when projected on the line join- 
ing MW to M31 is —165 ± 66, confirming the conventional 
views. On the other hand, the 3D inner motions inside the 
Local Group, if significant, are found to be surprisingly large 
(M31 subsystem has a relative 3D velocity of amplitude 399 
km s _1 with respect to the MW subsystem). It is interesting 
that the conclusion on the existence of transverse velocity 
of MW relative to M31 has been concluded previously by 
Peebles (1994) using its least action method (Peebles 1989). 
The point is discussed below in term of the convergence of 
the various luminous dipoles with the one observed in the 
Cosmic Microwave Background Radiation. 

The CMB temperature dipole, if interpreted as the sig- 
nature of our motion with respect to the rest frame of this 
radiation (CMB frame), gives for the Milky Way a pecu- 
liar velocity Vmw^cmb of 552 km s _1 pointing toward the 
galactic coordinates I — 266 and 6 = 29 (see Smoot et al. 
(1991, 1992), Kogut et al. (1993)). This MW motion in the 
CMB frame is traditionnally split into two components 

Vmw—cmb = Vmw-ilg + Vlg^cmb (5) 

where Vmw-lg is the velocity of the Milky Way relative 
to Local Group rest frame, which originates from the in- 
ternal non-linear dynamics governing the Local Group and 
Vlg^cmb is the peculiar motion of the Local Group as a 
whole in the CMB frame, created by large scale mass fluc- 
tuations present in the Universe. 

The latter can be infered in some way from the vari- 
ous luminous dipoles found in the literature (X-ray galaxies 
dipole by Miyaji & Boldt (1990): / = 313 and b = 38; Op- 
tical galaxies within 8000 km s _1 dipole by Hudson (1993): 
/ = 242 and b = 49 or I = 231 and b = 40; IRAS dipole, 
the least for shell (142.8-157.3 Mpc) by Plionis, Coles and 
Catelan (1993): I = 260.7 and b = 39.1; Abell/ACO clusters 
within 20 000 km s _1 by Branchini&Plionis (1996): I = 265 
and b = 16 if corrected from Virgocentric flow). Though 
there is no absolute convergence, it is noticed that conver- 
gence level is improved if deeper samples of extragalactic 
objects are involved. This fact indicates at least that X-ray, 
optical and IR observed objects have approximatively the 
same spatial distribution. 

Conversion of luminous dipole in terms of the motion 
Vlg^cmb of the Local Group in the CMB frame assumes 



the following hypotheses : 

- HI) The luminous objects trace the large scale mass 
density field. 

- H2) The linear approximation holds (in particular, 
peculiar velocity remains parallel to acceleration troughout 
the evolution of large scale structures). 

- H3) The sample of objects under consideration is at 
rest in the CMB frame. 

In what extent assumptions HI and H2 are satisfied is yet an 
open question, while results obtained from peculiar velocity 
analysis seem to challenge assumption H3 (Bulk flow of the 
shell of galaxies (within 3500-6500 km s _1 ) by Rubin et al. 

(1976) : t> B =950 km s _1 toward I = 308 and b = 25 in the 
CMB frame; Bulk flow of clusters of galaxies (up to 15 000 
km s _1 ) by Lauer&Postman (1994): i>.b=700 km s _1 toward 
I = 340 and b = 50 in the CMB frame, recently revisited by 
Graham (1996): t> B =738 km s _1 toward / = 330 and b = 45). 

On the other hand, the peculiar motion of the Milky 
Way in the Local Group rest frame Vmw^lg can be ob- 
tained from the analysis of the dynamics of the Local 
Group and its substructures. Substracting it to the observed 
Vmw^cmb thus gives a local estimate of Vlg-*cmb, which 
can be compared with its values extracted from dipoles anal- 
ysis, as mentioned above. 

Herein, a rough kinematical model is adopted, assigning 
to the Local Group the mean motion of its main substruc- 
tures (i.e. MW and M31 subsystems, equally weighted). It 
gives Vmw-,lg = (7 ± 303, 15 ± 155, -177 ± 144) or 178 
km s _1 toward I — 65 and b = —85. Finally, our local esti- 
mate of the LG peculiar velocity in the CMB frame yields 
Vlg^cmb = (-41 ± 303, -497 ± 155, 445 ± 144) or 668 km 
s _1 in amplitude pointing toward I = 265 and b = 42. This 
result can be directly compared with the well-known esti- 
mation by Yahil et al. (1977), based on the Solar system 
motion relative to the LG centroid : Vlg->cmb = 622 km 
s _1 toward I — 277 and b = 30 or (66, —535, 311) in galactic 
cartesian coordinates. 

Thus we have obtained the LG motion in CMB frame 
in an alternative way. Therefore it is remarkable that this 
result is in good agreement with the result of Yahil et al 

(1977) within 1 a level. The existence of some discrepancy 
between these two values, if significant, can be interpreted 
as follows. The LG centroid had been defined by specific 
choice of the main and satellite populations and with further 
search of the best-fit solution for the Solar system motion. 
In the present analysis we have found statistically significant 
indication of the bulk flow of the two main subsystems of 
LG, which can influence the definition of its centroid, and 
hence, the final result. 

The understanding of the cause of each discrepancy is 
the main problem to be solved. The LG substructure's 3D 
dynamics, as discussed above, could be essential for that 
problem. The fact of the existence of the bulk flow of the 
substructures should be crucial also while studing the past 
and future evolution of the Local Group. 

Moreover, bulk flows can be common properties of sub- 
groups of clusters of galaxies (Gurzadyan&Mazure, 1997) - 
galaxy associations, thus reflecting the role of merging and 
other basic trends in the formation mechanisms of the clus- 
ters. 
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APPENDIX A : THE BULK FLOW STATISTIC 

In this appendix, we use the maximum likelihood technique 
in order to derive the statistics of the 3 components of the 
bulk flow. Quantities, such as the variance of these estima- 
tors or the velocity dispersion estimate, are obtained by us- 
ing numerical simulations. 

We define the statistical model as follows. It is assumed, 
that the data sample consists of N independent objects, 
which follow the theoretical probability density dPth of Eq. 
(2): 

dPth = 4>(vi, v 2, V3) dvi dv2 dv3 x p(r, l,b) r 2 cos b dr dl db 

where 4>(vi , V2 , V3) is the 3D velocity distribution function 
(in galactic cartesian coordinates) and p(r, I, b) is the 3D 
spatial distribution function (in galactic coordinates). Here- 
after, we assume that the 3D velocity distribution function 



of the sample adopts the following form 

4>(v!,v 2 ,V3) = g(vi;v x ,a v ) x g(v 2 ;v y ,a v ) x g(v 3 ; v z , a v ) 

where g(x; xo, a x ) is a Gaussian of mean xo and of dispersion 
a x . The velocity field of the sample is thus split into a mean 
3D velocity vs = (v x ,v y ,v z ) (i.e. a bulk flow) plus a 3D 
random component, isotropic Maxwellian, centered on and 
of velocity dispersion cr v . 

For a given object, the observables are the galactic lon- 
gitude I and latitude b and the radial velocity v r given in 
the galactocentric frame (MWR frame) : 



v r — v.f = v-i cos I cos b + v 2 sin / cos b + V3 sin b 



(1) 



where f = (cos I cos b, sin I cos b, sin b) = (ri,f2,r 3 ) is the 
line-of-sight direction and v = (vi, 1)2, vz) the 3D velocity of 
the galaxy in the MWR frame. By successively integrating 
over the distance r and over 2 components of the 3D veloc- 
ity (say, of «i and V2), we express the observed probability 
density dPobs in terms of the observables : 

dPobs = g{v r ', v x cos I cos b + v y sin I cos b + v z sin b, a v ) dv r 



x rj(l, b) cosbdl db 



(2) 



where r/(l,b) is the projected angular distribution func- 
tion of the objects under consideration (i.e. rj(l, b) = 
J p(r,l,b) r 2 dr). Note, that the information on the dis- 
tances r of the galaxies is not used. For a data sample 
of N galaxies with measured {I ,b ,Vr}k=i,N , the efficient 
part of the natural logarithm of the likelihood function 
C = C{v x , v y ,v z ,a v ) reads thus : 



C 



In 



iV 



(V* -V x f \ -Vyr\ -Vzf\) ' 



(3) 



where f\ = cos l k cos b k , f\ = sinZ fc cosfe fc and f\ = sinfc* 
(see Eq. (|l])). Maximizing C with respect to v x , v y and v z 
gives the following set of equations 

(fl{v r - V x f! - V y T2 - V z f 3 )) = 
(f2{v r - V x fl - V y f 2 - V z f'j,)) = 
V z f 3 )) 



(4) 



{f 3 {Vr 



v x n - v y r 2 







where (.) denotes the average on the sample (for example 
(v r f 2 ) = UNJ2k=l v r^ = l/Ar£j V =i t v fe sinZ fe cos& fe ). This 
set of linear equations can be rewritten 



(•») 





v x 




(v r fl) 


M 


Vy 















where M is the 3x3 symmetric matrix 



M 



(fifi) (fif 2 ) (nf 3 ) 

{ri?2) (r2V2) (f2fs) 
(fif 3 ) (f 2 f 3 ) (f 3 f 3 ) 



(6) 



Finally, the unbiased statistics vs = (v x ,v y ,v z ) of the 3D 
bulk flow vs = (v x , v y , v z ) is obtained by inverting the 3x3 
matrix M : 



(7) 



v x 




(Vrh) 


Vy 


= rvr 1 


{v r r~2) 


. v z _ 
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This bulk flow statistic is robust. It does not depend on the 
distance of galaxies nor on the projected angular distribu- 
tion function rj(l, b) of the sample. However, the accuracy of 
the v x , v y and v z estimators depends on r/(l,b) and on the 
velocity dispersion a v of the dynamical system. We thus, 
compute the variance of these estimators by using numeri- 
cal simulations which are supposed to mimic the real data. 

The derivation of the velocity dispersion estimator a v 
is not straightforward. Maximizing the efficient part of the 
likelihood function £ of Eq. ([]) with respect to a v , we have 
the following equation 

at = ((v r - v x ri - v y f 2 - v z f 3 ) 2 ) (8) 

Unfortunatly, the velocity dispersion estimator <r* : 

a* 2 = ((v r — v x fi - Vyf 2 - v z f s ) 2 ) (9) 

obtained while replacing the bulk flow vs by its estimate 
v_b, is biased. The reasons of this bias are twofold. On one 
hand, the variance of the v x , v y and v z estimates leads to 
enhance the value of a* and thus, to overestimate the veloc- 
ity dispersion a v . On the other hand, because of the finite 
size of the sample, the random variables v x , v y , v z and the 
radial velocity v r of each objects are correlated. This feature 
contributes to underestimate a v when using the <rj statistic. 
As a matter of fact, the unbiased estimator of the velocity 
dispersion depends on the accuracy of the bulk flow estimate 
and thus on the velocity dispersion itself. In this paper, we 
have applied a computational iterative process on numerical 
simulations which furnishes an unbiased value for <j v . 

This paper has been produced using the Royal Astronomical 
Society/Blackwell Science M?gX style file. 
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